Pressure induced magnetic phase separation in Lag 75Cao.25Mn03 manganite. 
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The pressure dependence of the Curie temperature Tc(-P) in Lao.75Cao.25Mn03 was determined 
by neutron difTraction up to 8 GPa, and compared with the metallization temperature Tim{P) @]. 
The behavior of the two temperatures appears similar over the whole pressure range suggesting a key 
role of magnetic double exchange also in the pressure regime where the superexchange interaction 
is dominant. Coexistence of antiferromagnetic and ferromagnetic peaks at high pressure and low 
temperature indicates a phase separated regime which is well reproduced with a dynamical mean- 
field calculation for a simplified model. A new P-T phase diagram has been proposed on the basis 
of the whole set of experimental data. 

PACS numbers: 75.47.Lx, 62.50.-p, 64.75. Qr 



Perovskite-like manganites Lai_xCaa;Mn03 (LaCax) 
exhibit a variety of physical properties depending on the 
Ca concentration x. The strong correlation among mag- 
netic, electronic, orbital and transport properties of man- 
ganites makes these systems particulary sensitive to ex- 
ternal perturbations, such as temperature variation, ap- 
plication of magnetic field or high pressure (HP) [ll, y]- 
The most popular phenomenon is the Colossal Magneto- 
Resistance (CMR) 01 observed over the 0.2 < x < 0.5 
range. Here, the system shows a transition from a high 
temperature paramagnetic (PM) insulating phase to a 
low temperature ferromagnetic (FM) metallic one and 
its properties can be described by two competing mech- 
anisms: double-exchange (DE) [J] mechanism and Jahn- 
Teller (JT) effect 0]. The DE model, which couples mag- 
netic and electronic degrees of freedom, qualitatively ex- 
plains the observed transition and the contiguity between 
the insulator to metal transition temperature T/j\/ and 
the Curie temperature Tc- More realistic modeling is 
only obtained by introducing the JT mechanism with its 
charge localizing effect [sj. 

HP techniques were successfully employed to investi- 
gate the role of JT structural distortion on manganite 
properties. The idea is that lattice compression reduces 
JT distortion 0,01, increases hopping integral t, thus fa- 
voring the onset of metallic phase |8| . Several manganites 
show indeed an increase of Tj m on applying pressure up 
to 2-3 GPa Nevertheless, experiments at higher 

pressures cast doubts on the above scheme 113 • Sev- 
eral HP experiments were carried out on LaCax mangan- 
ites. In particular, in LaCaO.25 Tim{P) was found to in- 
crease linearly only up to P~ 3 GPa, to bend down above 
6 GPa and to approach an asymptotic value close to room 
temperature [1, [T^l- An antiferromagnetic (AFM) su- 
perexchange (SE) interaction between the Mn magnetic 
moments [14| . which are formed by localized t2g elec- 



trons, was proposed to be responsible for the anomalous 
T/A/(^') evolution at HP [1^. This peculiar pressure be- 
havior can be ascribed to a competition between the DE 
and the SE interaction, which, is responsible for AFM 
order observed away from the CMR concentration re- 
gion at ambient pressure. Since the strength of DE is 
proportional to the hopping integral t and SE is pro- 
portional to i^, on increasing pressure (i.e. t) the role 
of SE versus DE is continuously enhanced and it is ex- 
pected to become the dominant interaction in the HP 
regime [l5| . The pressure dependence of the magnetic 
transition was recently investigated by neutron diffrac- 
tion in LaCaO.33 ^ and LaCaO.25 ^ over the - 4 
GPa range and an almost linear P-dependence of Tc was 
observed. In LaCaO.25 rather close values of Tjm and Tc 
[l7| were found over the whole P range. As a confirm of 
the role of SE, evidence of AFM order coexisting within 
the FM phase was also obtained at low temperature and 
HP [l^, [l3 • The suppression of the long range FM or- 
der at HP (around 23 GPa) in LaCaO.25 was deduced 
by X-ray magnetic dichroism data [l8| . The remarkable 
compression of MnOg octahedra along the b axis and the 
anomalous change in the lattice strain at around 23 GPa 
together with the monotonic decrease of the FM mo- 
ment with pressure, suggested indeed a transformation 
from an orbitally disordered FM phase to an orbitally or- 
dered AFM phase [l^ ■ The coexistence of the two mag- 
netic phases in the pressure range between 2 and 23 GPa 
was conjectured [131. It is worth to note that although 
the tendency toward mixed phases, whether structural or 
magnetic, was recently observed in several manganites at 
ambient [l^ and high pressure @, [20l - l22l | , the theoretical 
investigation of this phenomenon was rather limited. 

Although LaCaO.25 is the most extensively investi- 
gated manganite system under pressure, no experimen- 
tal data are available on Tc{P) over the P range where 
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Tim is no more linear and the SE mechanism is dominant 
over the DE. In the present paper we extend the neutron 
diffraction measurements on LaCaO.25 up to 8 GPa in 
order to clarify the role and the relative strength of DE 
and SE interactions and the coupling between and 
Tim hi HP regime. The tendency towards phase separa- 
tion is modeled by using the simple theoretical approach 
successfully employed in Ref. |15| . 

The neutron diffraction experiment was carried out at 
the LLB (Saclay, FR) using the G6 — 1 two-axis pow- 
der difFractometcr with a selected incident neutron wave- 
length A = 4.74 A. The diffraction patterns were col- 
lected over the 1.5-300 K temperature range along sev- 
eral isobaric paths. Details about the cell are reported in 
Ref. Sapphire anvil were used at 2.1, 3.9, and 6.0 

GPa and, diamond anvils were employed at 8 GPa. The 
sample volume ranged from several millimeters (sapphire 
anvils) to 0.01 mm'^ (diamond anvils). The sample, pre- 
pared by a solid reaction method, was finely milled and a 
mixture of sample and NaCl salt (1:2 volume proportion) 
was loaded in the gasket hole together with small ruby 
spheres for pressure calibration. 

Ambient pressure data were collected over a wide d- 
space range 2.5 A < d < 54 A (20 = 5 - 140°) using two 
detector positions. Diffraction patterns at high and low 
temperature are shown in the inset of Fig[lja). At low 
temperature, the strong increase of the peak intensities at 
d=3.87 A ((1 l)/(0 2 0)) and d=2.74 A ((2 OYO 2)(1 
2 1)) can be ascribed to the FM contribution [131 (inset 
of FiglTJa)). The diffraction pattern was thus refined us- 
ing the Pnma orthorhombic structure with an associated 
FM order. The lattice parameters (a=5.47 A, 6=7.72 A, 
c=5.50 A) and the magnetic moment fj, = 3.57/iB ob- 
tained at ambient pressure by Rietveld refinement are in 
good agreement with the results of Ref. ^J] . 

In order to reduce the acquisition time at HP, we fo- 
cused on the main peak at d=3.87A and diffraction pat- 
terns were collected exploiting only one detector position 
{29 = 12 - 91° and 3.3 A < d < 22 A). Data collected at 
6 GPa are shown in FiglUb) where new magnetic reflec- 
tions are found at d=7.38 A (0 1 0) and d=3.47 A (1 1 1). 
These peaks are unambiguously observed at P= 3.9 GPa 
and upwards, whereas the data collected at P= 2.1 does 
not allow any definitive conclusion. Nevertheless, a new 
pressure induced peak (d=7.47 A at 1.5 GPa) was previ- 
ously reported in Ref. [ij] and ascribed to an AFM order 
with A-type structure and propagation vector k — (000). 
This suggests that the AFM peak is probably too weak 
compared with the signal/noise ratio to be observed in 
the patten collected at P= 2.1 GPa. 

Looking at the present data, a discrepancy is observed 
between the observed AFM peak position and the ex- 
pected rf-spacing calculated on the basis of the lattice 
parameters. For example, considering the lattice param- 
eters ((a=5.4lA, 6=7.66A c= 5.45A).) obtained at 6 
GPa, the (0 1 0) AFM peak is expected to be found at 
7.66 A, whereas it is centered at d=7.38 A (see FiglT](b)). 
At this point, it is worth noticing that both Kozenklo and 
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FIG. 1: (a): Diffraction pattern and Rietveld refinement 
at P=0 and T=1.5 K. Inset: comparison between the P=0 
diffraction patterns collected at T=240 and 1.5 K. (b) Diffrac- 
tion patterns at 6 GPa collected at T=280 and 1.5 K. Peaks 
at d = 3.47 and d = 7.38 A are ascribed to AFM order. Cal- 
culated positions of the Bragg reflections are also shown: the 
two upper sets correspond to the nuclear and magnetic phase 
in agreement with Ref. [l3, , the two lower sets are asso- 
ciated to the second structural and magnetic phase. Inset: 
FWHM of AFM and nuclear peaks versus pressure at T=1.5 
K. 



Ding observed an uniaxial anisotropic compression of the 
MnOg octahedra along the b axis of the orthorhombic 
structure, which is consistent with this shift of the AFM 
peak position [13, [13 • These findings suggest that the 
two magnetic phases, the AFM and the FM, are not as- 
sociated to the same nuclear phase. Two possible expla- 
nations can be proposed: i) the AFM order is not com- 
mensurate with the nuclear phase, as already reported 
for some manganite compounds [23I; ii) the AFM order 
is associated to a nuclear phase with the same Pnma 
symmetry but different values of lattice parameters. 

In support of the ii) hypothesis, the positions of nuclear 
and magnetic reflections were obtained above 4 GPa us- 
ing Fullprof software (see Figlljb)) [2^. For example at 
6 GPa, two Pnma nuclear structures and two magnetic 
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FIG. 2: Normalized magnetic moments versus temperature 
(solid lines are guides to the eye). Inset: temperature depen- 
dence of the FM peak h,k,l=l,0,l 



FIG. 3: Temperature variation of Tc and Tim 0] versus 
pressure. Inset: tentative P-T phase diagram. The onset of 
the AFM peak at 2 GPa refers to Kozenklo et al. work [13] 



orders were considered: the first one associated to FM 
order with the set of lattice parameters reported in Ref. 
pTl fisj ) and the second one associated to an AFM or- 
der with a=5.58 A, 6=7.38 A and c=5.49 A. The AFM 
unit cell was described by four Mn atoms: Mui =(00^), 
Mn2 = (iOO) Mna = (Qii) Mn4 = (^^0). The magnetic 
structure is a combination of the and Fy representa- 
tions. Bragg reflections calculated for the AFM structure 

(A-type) using the R^; = (-1- H ) and Ry = (+ H ) 

sequence of magnetic moments and propagation vector 
k = (000), were well consistent with the experiment as 
shown in Figlljb). The magnetic moments lay on the 
ah plane and they are placed in a way that is similar to 
that one observed in LaMnOs at high pressure [13]. The 
fractional volumes of the FM and the AFM phases were 
obtained for the neutron pattern collected at T= 5K and 
were found to pass from 87.9 % (FM) and 12.1 % (AFM) 
at 4 GPa to 76.3 % (FM) and 23 % (AFM) at 8 GPa. The 
FWHM of the FM (101)/(020), and the AFM peak (010), 
show opposite behaviors on increasing the pressure: the 
former increases while the latter decreases (see the inset 
of Fig(IJb)). This behavior is thus consistent with the 
hypothesis ii), that is lattice compression simultaneously 
induces both AFM order and a separation between two 
isostructural nuclear phases. 

The emerging scenario is consistent with a founded 
physical picture. If AFM order is driven by SE inter- 
action, Mn orbitals are expected to order in a way that 
reduces h axis. The connection between lattice and or- 
bital degrees of freedom as investigated in several or- 
thorhombic manganites [2ll [28l . [29| . suggests that the 
development of orbital ordering results in a contraction 
of h parameter, and if h/y/2 < c < a the occurrence 
of an orbital ordering can be conjectured [1^, H^. On 
the other hand, if c > a ~ b/\/2, orbital disorder is ex- 
pected. It appears thus reasonable to associate FM and 
AFM order to orthorhombic phases with large and small 



b parameters: 7.66 A and 7.38 A respectively. 

The FM to PM transition can be identified looking 
at the temperature dependence of the integrated inten- 
sity of the main peak (d=3.87A) (inset in Fig. [5]). The 
results of the refinements do not considerably change in- 
cluding or excluding the Debye- Waller factor suggesting 
that the temperature dependence of the nuclear intensi- 
ties can be neglected. Therefore, the Debye- Waller factor 
was neglected over the scattering vector range relevant 
to magnetic peaks. After proper background correction, 
the integrated peak intensities, subtracted by the nuclear 
contribution, can be reasonably considered proportional 
to the effective magnetic moments M. The temperature 
dependence of the normalized magnetic moments M/Mp 
(where Mq is the effective magnetic moment at the lowest 
T) is shown in Fig [5] at different pressures. 

A rapid increase of Tc on applying pressure is evident 
over the low pressure range (0-4 GPa) , whereas a strongly 
reduced pressure dependence is observed above 4 GPa. 
We notice that the observed broadening of the transition 
can be likely ascribed to strain effects due to pressure 
gradients. The data were analyzed using a modified Bril- 
louin function to describe the temperature dependence 
of the magnetization, according to the Weiss molecular 
field theory [30l |. To avoid the intrinsic ambiguity of the 
definition of both Tc and Tjm, we used the tempera- 
ture variations AT from ambient pressure Tc(P=0) and 
T/m(P=0) to compare the pressure dependence of Tc(P) 
and T/jv/(P) [I]- As shown in Fig. [31 the pressure de- 
pendence of AT is basically the same for both Tc and 
T/j\/: a linear increase at low pressure followed by a sat- 
uration towards an asymptotic value, above 6 GPa. This 
result further confirms the onset of a pressure driven lo- 
calizing mechanism that, becoming competitive with DE 
delocalizing mechanism at high pressure, strongly affects 
both the transition temperatures Tc and Tjm- More- 
over, the onset of an AFM state at low temperature co- 
existing with the FM order confirms the SE-AFM nature 
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of the localizing interaction. This suggests a phase sep- 
aration scenario at low temperature, as depicted in the 
inset of Fig. |3] where a tentative P-T phase diagram 
is shown. We finally notice that AT c{P) slightly de- 
viates from AT/m(-P) only at very high pressure. This 
means that decreasing temperature at a pressure above 
^ 4 GPa, LaCaO.25 enters a metallic phase before achiev- 
ing the FM state. This features might be consistent with 
a percolation scenario among metallic clusters in a PM 
insulating bulk, as previously suggested by HP optical 
measurements [20|. 

In order to understand theoretically if the pressure- 
induced increase of AFM-SE can, in turn, induce a phase 
separation, we considered a simplified model for mangan- 
ites similar to that used in Refs. [TBI. [3T|. 

(•y)cr i<7 

—Jh ^ ■ Si + 

i 

+Ji^Si ■ Sj - (rii) {ai + al) +ujoa'^a.{l) 

(ij) i 

A single band of itinerant electrons with nearest- 
neighbor hopping t replaces the two eg bands, and a Hol- 
stein coupling to a local Einstein phonon of frequency 
luq mimics the JT interaction. The DE physics is in- 
troduced through an FM coupling between the itinerant 
band and localized spins Si associated to the t2g orbitals. 
The localized spin are, in turn, antiferromagnetically cou- 
pled with a coupling constant which is small at ambient 
pressure and grows with increasing pressure ad discussed 
in Ref. [l^. We work in the grandcanonical ensemble, 
where the chemical potential fi determines the particle 
density in the Cg bands. The density n = (^ig- ni_a-) , 
where rii.o- = cl gCi,^ is related to the Ca doping by the 
relation x = 1 — n 

The study of the thermodynamic instability associated 
to phase separation requires to work in the thermody- 
namic limit, which can be addressed using dynamical 
mean- field theory (DMFT)!!^. We studied the Holstein- 
DE model (H]) at T=0 using DMFT and exact diagonal- 
ization to solve the numerical part of DMFT. We used our 
implementation of DMFT, as discussed in Ref. [s^, with 
exact diagonalization using 8 levels in the discretized 
bath. 

Since the aim of this calculation is to understand if, on 
general grounds, phase separation can be expected from 
the competition between AFM and FM, we introduced 
some safe approximations. In particular we treated the 
localize t2g spins as classical variables neglecting quan- 
tum fluctuations. Moreover, in light of the experimental 
suggestions, only the two configurations FM and com- 
mensurate AFM ordering were considered, neglecting in- 
termediate orientations and more exotic spin patterns. 
We computed the total energy for these two phases and 
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FIG. 4: P-x phase diagram for the Holstein-DE model with 
parameters from Ref. [l^. A pressure dependence line sepa- 
rates a phase separation region from a homogeneous ferromag- 
netic state. Phase separation is favored by increasing pres- 
sure. Insets: the x-dependence of the Free Energy normalized 
to the bandwidth at P=0 and 6 GPa. Dashed red (blue) line 
represents the x-dependence for the AFM (FM) phase. Solid 
black line is a Maxwell construction whose boundaries deter- 
mine the doping range for which phase separation takes place 
at the chosen pressure. 



for the nonmagnetic phase as a function of the Ca con- 
centration. The effect of pressure has been included in 
the calculations adopting the parameters of Ref. (l5l |. 
with the SE coupling scaling as the square of the hopping 
parameter. Even though our model is too simplified to 
attempt a quantitative comparison with the experimental 
phase diagram, it is expected to properly reproduce the 
competition between magnetic phase which ultimately 
leads to phase separation. 

Phase separation is associated to a negative charge 
compressibility k = dn/d^, which corresponds to a neg- 
ative curvature of the energy as a function of density 
(or equivalently on doping). In order to detect phase 
separation we need to solve the system for various val- 
ues of Ca concentration and analyze the curvature of the 
ground state energy around the experimentally relevant 
case x = 0.25. In our case each of the independent solu- 
tions (FM, AFM, non magnetic) has the correct positive 
curvature for every value of x (see curves in the two inset 
of Fig. |4|). However, when two curves cross as a function 
of doping and the system undergoes a first-order transi- 
tion from one phase to the other, the energy has a cusp, 
and. more importantly, one can lower the energy by di- 
viding the system into a fraction of AFM and a fraction 
of FM. This is realized through a Maxwell construction 
(solid black line in the two panels) and physically corre- 
sponds precisely to phase separation. 

In FigUlwe show the pressure vs Ca doping phase di- 
agram obtained with DMFT. 

At zero pressure we find phase separation between 
AFM and FM only close to a; = [SJI- Our simpli- 
fied model shows indeed phase separation below x ~ 
0.16, while larger doping correspond to homogeneous FM 
phases. Increasing the pressure increases the phase sep- 
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arated region, pushing its boundary to larger doping. At 
P ~ 5 GPa the phase separation region reaches x = 0.25. 
Further increasing pressure at fixed x enhances the size of 
the AFM fraction of the system. Despite the simphfica- 
tions introduced in our model, the characteristic pressure 
leading to phase separation is reproduced accurately by 
the theoretical calculations, confirming that the enhance- 
ment of the AFM superexchangc can indeed account for 
the observed pressure-driven phase separation. 

In summary, the pressure dependence of the magnetic 
temperature Tc(P) was obtained up to 8 GPa and 
found to be rather close to Tim{P), suggesting that 
DE model still holds at HP. The reinforced role of SE 
localizing mechanism, competing with delocalizing DE, 
was confirmed by Tc{P) and by the onset of an AFM 
state at low temperature. The novel antiferromagnetic 
structure is not related to the crystalline structure 
associated to the FM phase and observed at ambient 
pressure. The data are consistent with a picture in 



which lattice compression simultaneously induces both 
AFM order and a separation between two isostructural 
nuclear phases. The simultaneous appearance of FM 
and AFM peaks depicts a phase separation scenario 
promoted by pressure. This tendency towards a phase 
separation is properly reproduced, using a simplified 
theoretical model indicating that formation of domains 
is promoted by pressure. 



Acknowledgments 

M.C. is funded by the European Research Council 
under FP7/ERC Starting Independent Research Grant 
"SUPERBAD" (Grant Agreement n. 240524) and MIUR 
PRIN 2007 Prot. 2007FW3MJX003. 



[1] E. Dagotto et al. Physics Reports 344,1-153 (2001) [20 
[2] S.-W. Cheong and H. Y. Hwang, in Colossal Magnetore- [21 
sistance Oxides, edited by Y. Tokura, Monographs in [22 
Condensed Matter Science (Gordon and Breach, Read- [23 
ing, U.K., 2000). [24 
[3] S. Jin et al, Science 264, 413 (1994) [25 
[4] C. Zener, Phys. Rev. 82, 403 (1951) [26^ 
[5] A. J. Millis, Nature 392, 147 (1998) (London); A. J. Millis [27 

et al, Phys. Rev. Lett. 74, 5144 (1995). 
[6] I. Loa et al, Phys. Rev. Lett. 87,125501 (2001) [28 
[7] A. Congeduti et al. Phis. Rev. Lett. 86, 1251 (2001) [29 
[8] P. Postorino et al, Phys. Rev. Lett. 91, 175501 (2003) [30 
[9] H. Y. Hwang et al, Phys. Rev. B 52, 15046 (1995) [31 
[10] V. Laukhin et al, Phys. Rev. B 56 R10009(1997) 
[11] C. Cui et al, Phys. Rev. B 68 214417 (2003) [32 
[12] C. Cui et al, Phys. Rev. B 70 094409 (2004) [33 
[13] A. Congeduti et al, Phys. Rev. B 63, 184410 (2001). 
[14] D. Feinberg et al, Phys. Rev. B 57 R5583 (1998); M. [34 
Capone et al, Europhys. J. B 17,103 (2000) [35 
[15] A Sacchetti, et al. New J. Phys. 8, 3 (2006) 
[16] D. P Kozlenko et al, J. Magn. Magn. Mater. 258, 290 
(2003) 

[17] D. P. Kozlenko et al JETP Letters 82, 447 (2005) 
[18] Y. Ding et al, Phys. Rev. Lett. 102, 237201 (2009) 
[19] E. Dagotto, New J. Phys. 7, 67 (2005) 



A. Sacchetti et al, Phys. Rev. Lett. 96, 035503 (2006) 
L. Malavasi et al, App. Phys. Lett. 94, 061907 (2009) 
M. Baldini et al, Phys. Rev. B 80, 045123 (2009) 
I. N. Goncharenko, High Press. Res. 24, 193 (2004) 
P. G. RadaeUi et al, Phys Rev. B 56, 8265 (1997) 
S. Lan dsgesel l et al, Phys. Rev. B 80, 014412 (2009) 
|http: / /www.ilTeu/sites /fullprof/ 1 

L. Pinsard-Gaudart et al, Phys. Rev. B 64, 064426 
(2001) 

P. M. Woodward et al, Chem. Mater. 10, 3652 (1998) 
L. Malavasi et al, J. Mater. Chem. 10,1304 (2010) 
A. B. Beznosov et al, Phys. Rev. B 68, 054109 (2003). 
M. Capone and S. Ciuchi, Phys. Rev. B, 65, 104409 
(2002). 

A. Georges et al. Rev. Mod. Phys. 68, 13 (1996) 

M. Capone, L. de' Medici, A. Georges, Phys. Rev. B 76, 

245116 (2007). 

R. S. Fishman et al, New J. Phys. 8, 116 (2006). 
Details of the calculations and a comparison with more 
realistic models of the manganites (inclusion for orbital 
degeneracy of the eg electrons and proper JT interaction, 
quantum description of t2g spins) will be presented in a 
forthcoming work. 



